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ABSTRACT 

We revisit the Polish doughnut model of accretion disks providing a comprehensive 
analytical description of the Polish doughnut structure. We describe a perfect fluid cir- 
cularly orbiting around a Schwarzschild black hole, source of the gravitational field, by 
the effective potential approach for the exact gravitational and centrifugal effects. This 
analysis leads to a detailed, analytical description of the accretion disk, its toroidal 
surface, the thickness, the distance from the source. We determine the variation of 
these features with the effective potential and the fluid angular momentum. Many an- 
alytical formulas are given. In particular it turns out that the distance from the source 
of the inner surface of the torus increases with increasing fluid angular momentum but 
decreases with increasing energy function defined as the value of the effective poten- 
tial for that momentum. The location of torus maximum thickness moves towards the 
external regions of the surface with increasing angular momentum, until it reaches a 
maximum an then decreases. Assuming a polytropic equation of state we investigate 
some specific cases. 
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1 INTRODUCTION 

Accretion disks are one of the most intriguing issues in high energy Astrophysics. They enter into very different contexts: 
from the proto-planetary disks to Gamma Ray Bursts (GRB), from X-ray binaries to the Active Galactic Nuclei (AGN). 
Indeed, many aspects of disk structure, its dynamics, the formation of "jets", the exact mechanism behind the accretion, its 
equilibrium, the confinement and stability under perturbations are still uncertain. 

Analytic and semi-analytic models for accretion onto a compact object are generally stationary and axially symmetric. 
Thus, all physical quantities depend only on the radial distance from the center r, and the vertical distance from the equatorial 
symmetry plane z. Thick accretion disk (z/r ^> 1) models are basically designed according to the following assumptions: first, 
a circular motion dynamic is prescribed. Second, the matter being accreted is described by a perfect fluid energy-momentum 
tensor. This requirement lies on the assumption that the time scale of dynamic processes, involving the pressure, centrifugal 
and gravitational forces, are smaller, and much smaller in the case of vertically thin disks, than thermal ones, which in turn 
are smaller or much smaller than the viscous time scale. This implies that, at variance with many thin disk models, dissipative 
effects like viscosity or resistivity are neglected, and accretion is a consequence of the strong gravitational field of the attractor. 
Indeed this is a great advantage, since angular momentum transport in the fluid is perhaps one of the most controversial 
aspects in thin disk models. 

In this work we focus on the Polish doughnut model. The Polish doughnut is a fully relativistic model of thick accretion 
disk with a toroidal shape, and is an example of opaque (large optical depth) and super-Eddington (hight matter accretion 
rates) disk model. During the evolution of dynamic processes, the functional form of the angular momentum and entropy 
distribution depends on the initial conditions of the system and on the details of the dissipative processes. Paczyriski realized 
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that it is physically reasonable to assume ad hoc distributions. The Polish doughnut is characterized by a constant angular 
momentum ( |Abramowicz||2008[ ) . 



The development of this model was drawn up by Paczyhski and his collaborators in a series of works (Paczyhski 



Paczyhski & Wiita 1980; Kozlo wski et al. 1978; Ab ramowicz et al.|197 8; Jaroszynski et al. 1980; Abr amowicz et al.|1980l 



1980 



1997 



Fishbonefe Moncrief|1976| ). Using a perfect fluid energy- momentum tensor, Abramowicz et aT] ( |1978| wrote the equations of the 



hydrodynamics for this model. Jaroszynski et al. (1980) discussed the case of a non barotropic fluid and showed an important 
result concerning the pressure of the rotating fluid: for a perfect fluid matter circularly rotating around a Schwarzschild black 
hole, the shapes and location of the equipressure surfaces follow directly from the assumed angular momentum distribution 
(Jaroszynski et al. 1980). More recently, |Lei et aL (2009) assumed an angular momentum distribution in the form that 
depends on three constant parameters, and different configurations have been studied. A significant result concerning this 
model, known as the von Zeipel condition, has been extensively investigated in |Kozlowski et al." ( 1978 ) and Jar oszynski et| 



al. (1980): the constant pressure surfaces coincide with surfaces of constant density if and only if the surfaces at constant 



angular momentum coincide with the surfaces at constant relativistic velocity. More generally we can say that in the static 
spacetimes the family of von Zeipel's surfaces does not depend on the particular rotation law of the fluid, in the sense that it 
does not depend on nothing but the background spacetime. An accurate study of the von Zeipel surfaces has been performed 
in |Abramowicz1 ( |1971t and |Chakrabarti| fl!99l| |1990 ) . 

Paczyhski realized from the study of Roche lobe in the accretion disks of the binary systems that the black hole Roche 
lobe overflow must induce the dynamical mass loss from the disk, thus the accretion (Boyer 1965; Frank et al. 2002). The 
accretion occurs at the point of cusp of equipotential surfaces. This process is realized by the relativistic Roche lobe overflow. 
This clearly is an explanation for the accretion that does not involve other factors (as the dissipative ones) than the strong 
gravitational field of the attractor. However Abramowicz (1981) showed that it constitutes also an important stabilizing 
mechanism against the thermal and viscous instabilities locally, and against the so called Papaloizou and Pringle instability 
globally flBlaes|1987| ). 

The general relativistic effects on matter dynamics close to a Schwarzschild black hole have been modeled in an approx- 
imate pseudo-Newtonian theory by Paczyhski and Wiita introducing a properly chosen non-exact gravitational potential, 
known as Paczyhski- Wiita (P-W) potential (Paczyhski & Wiita 1980). This potential simulates the relativistic effects of 
the gravitational field acting on the fluid in the disk, in the Schwarzschild spacetime. This is not the exact expression of 
the effective potential for gravitational and centrifugal effects, and yet the P-W potential cannot properly be considered a 
newtonian approximation, that is valid in the limit of weak gravitational fields. The P-W potential differs from the exact 
relativistic one by a constant in its gravitational part and it is a Newtonian way to write some of the general relativistic effects 
characterizing the thick disks. With this approximation, the radius of a marginally bound orbit, the last stable circular orbit 
radius of the Schwarzschild spacetime and the form of the Keplerian angular momentum have been correctively reproduced. 
A step-by-step derivation of the P-W potential and a detailed discussion of its main features can be found in |Abramowicz| 
(2009). The agreement between model predictions and simulations of accretion flows has been verified and found an excellent 
outcome, in e.g. Igumenshche & Abramowicz (2000) and Shafee et al. (2008). The equipressure surfaces for a Schwarzschild 



black hole have been compared with global magnetohydrodynamic numerical simulations in Fragile et al. (2007), (see also De 



Villiers fe Hawley||2002l |Hawley|[l987l p90| [1991] |Hawley et al.||1984|). Recently, the study of the Polish doughnut model 



has been developed for different attractors (see Stuchhk et al.||2009| | Stuchhk &; Slany||2006| Stuchhk &; Kovaf||2008| for the 
Schwarzschild-de Sitter and Kerr-de Sitter spacetimes). 

In the present work we face the study of the Polish doughnut model in the Schwarzschild background using the approach 
of the general relativistic effective potential in its exact form. Gravitational and centrifugal forces carried out in the effective 
potential and the pressure force operate on a perfect fluid of the disk. When the latter vanishes, the hydrodynamics of the 
fluid describes a geodesic disk whose equations are formally resembling those of motion of a test particle orbiting in the same 
background. We take advantage of this formal analogy using the familiar and well known results on the dynamics of the test 
particles to get a comparison between the Polish doughnut, which is supported by the pressure, and the geodetic disk. In this 
way we can evaluate the right weight of the pressure effects on the dynamics of fluid and the shape of the torus. Furthermore, 
we draw a complete and analytic description of the toroidal surface of the disk, including the analysis of its extension in 
space, the distance from the center attractor, its thickness etc, and understand how these features are modified by changing 
the angular momentum of the fluid and the effective potential. In particular we find that the distance from the source of the 
inner surface of the torus increases with increasing fluid angular momentum and decreases with increasing energy function 
defined as the value of the effective potential for that momentum. 

In Section [2] we introduce the Polish doughnut model writing the equations of the ideal hydrodynamics for a fluid 
circularly orbiting in the background of the Schwarzschild spacetime. In Section [3] we detail the fluid pressure gradients along 
the radial direction, and in Section [4] along the polar angular direction to determine the regions of maximum and minimum 
pressure in the disk, the regions of increasing pressure and the isobar surfaces. In Section [5] we trace the profile of the toroidal 
disk in the Polish model by introducing and studying in detail the Boyer potential for the barotropic fluid. In Section [6] we 
investigate the case of poly tropic equation of state. In Section [7] we analyze the proper fluid velocity, finding the regions of the 
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disks of maximum and minimum velocity. In Section [8] we discuss the relativistic angular velocity examining the properties 
of the von Zeipel surfaces. Conclusions follow. 



2 THE POLISH DOUGHNUT MODEL 

Consider a one-species particle perfect fluid (simple fluid), where 

T ab = (p + p)U a U b + Pgab (1) 

is the fluid energy momentum tensor, p and p are the total energy density and pressure, respectively, as measured by an 
observer moving with the fluid, and g a b the metric tensor. The time-like flow vector field U denotes the fluid four-velocit}|^] 
The motion of the fluid is described by the continuity equation: 

U a V aP +{p + p)V a U a = 0, (2) 

and the Euler equation: 

{p + p)U a V a U c + h bc V b p = 0, (3) 



where h a b = g a b + UJJb (Mi sner et al.|1973 ) 



Neglecting the fluid back reaction, we consider the fluid motion in the Schwarzschild spacetime background: 



-e 



/(r) dt 2 + e~ Hr) (r)dr 2 + r 2 (di9 2 + sin 2 tidtp 2 ) , (4) 
written in standard spherical coordinates, where = (1 — 2M/r). We define: 

A = U r , E = tf', $ = 17*, 9 = U" & , (5) 
and we introduce the set of variables {E, V sc , L, T} by the following relations: 

A = y/E* - Vh , 2=4' ®= 2 L 2 Q > e =^> (6) 
where 



V.c = <le»M I 1 + (7) 



r 2 sin 2 $ 

is the effective potential (Mis ner et al.|1973 ). In fact, from Eq. [6] we obtain: 



r 2 = (E 2 - Vl) (8) 

(the dot represents differentiation with respect to the proper time): Eq. [8] describes the motion inside the effective potential 
V sc , defined as the energy at which the (radial) kinetic energy of the fluid element vanishes. 

Eq.[8]and the definitions in Eqs. |5|6(7| are formally the same as for the test particle motion in the Schwarzschild spacetime. 
Obviously, for the particle motion, U a in Eq. [5] is the test particle four-velocity, and (E, L) in Eq. [6] are two constants of 
motion, the particle energy and angular momentum per unit of mass as seen by infinity, respectively. 

Using this similarity, we can make a one-to-one comparison of the motion of the fluid, under the action of the pressure 
forces balanced by the effective potential, with the test particles dynamics regulated by the gravitational and centrifugal forces 
as described by the effective potential (Wald 1984; Misner et al.||1973| Pugliese et al.|[2011a|b|c ). The comparison with the 



case of dust disk through the effective potential enables us to evaluate the relationship between the contribution of pressure 
and the gravitational and centrifugal effects to the dynamics of the system, especially in relation to the angular momentum 
of the fluid in rotation and the disk shape. 

We consider the case of a fluid circular configuration, defined by the constraints A = (i.e. V sc — E), restricted to a 
fixed plane sin^ = a/0. No motion is assumed in the $ angular direction, which means G = 0. For the symmetries of the 
problem, we always assume d t Q = and c^Q = 0, being Q a generic tensor of the spacetime (we can refer to this assumption 
as the condition of ideal hidrodynamics of equilibrium). 

Within our assumptions (A = 0, G = 0, dtp = d^p = 0), from the Euler equation [3] we derive the expressions for the radial 
pressure gradient G r and the angular pressure gradient G$: 

1 The fluid four-velocity satisfy U a U a = — 1. We adopt the geometrical units c = 1 = G and the (—,+,+,+) signature. The radius 
r has unit of mass [M], and the angular momentum units of [M] 2 , the velocities [{/*] = [U r ] = 1 and [U^] = [U^] = [M]~ 1 with 
[U^ /U*] = [M] _1 and [U^/Ut] = [M]. For the seek of convenience, we always consider the dimensionless energy and effective potential 
[V sc ] = 1 and an angular momentum per unit of mass [L]/[M] = [M]. 
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and 

G , = ^ = Wl^r 2 * 2 , (10) 
P + P 



3 THE RADIAL PRESSURE GRADIENT G 



R 



The first part of the present work is dedicated to the study of the fluid angular momentum, in particular, we are interested 
especially in the comparison between the geodetic disk case, and the case of a fluid subjected to a non-zero pressure, for this 
purpose in this and following Sections we will study the radial and angular pressure gradient. This study allows us to evaluate 
the pressure contribution to the disk dynamics along the orbital radius and the plans on which the accretion disk stretches. 
This analysis introduces the second part of the work in which we finally trace the profile of the disk obtained from the analysis 
of the constant pressure surfaces. 
Equation [5] can be written as: 



(12) 



p + p 2 \ dr 

in terms of the partial derivative of V sc computed keeping L = constant. Eq. [II] yields to: 

'dV? r ' 



d r p — for 



and 
lim G r 



0, 



dr 



lim G r 

^2M 



o, 



(13) 



Assuming p > and p > 0, from Eq. 11 it follows that the pressure increases (decreases) with the orbital radius r as V sc 
decreases (increases), and that the critical points of p (as a function of r) are the same as of V sc as a function of r at L 
constant. Thus, solving Eq. [l2]for the unknown L we find that these critical points are for: 



a 2 Mr 2 



(14) 

3M is the last circular orbit radius for a test particle in 
is therefore satisfied only in the range r > n co , for fixed a. In r = n co it is 
G r = — 1/(3M). The angular momentum L K describes the isobar fluid configurations: where the fluid is characterized by 
L = L K , the pressure p is constant and the Euler equation [3] reduces to U a V a U h — 0, describing the motion of a pressure- free 
fluid (dust). The curves L = L K represent the critical points of the pressure p. From Eq. 14 it follows, according to the physics 



-3M) " 

This function is defined in the range r > n co , where n 
the Schwarzschild spacetim^] Equation 



12 



exists in the range [2M,n co ], where d r p < (pressure always decreasing) ( Kozlowski et al.||1978 
Paczyhski||1980| |Abramowicz]|2009| |Lei et al.|2009 



of the free (test) particle (and dust defined by p = 0) ( |Misner et al.|[l973 Pugliese et al.||2011a|b|c ) that no critical point 



Abramowicz et al. 



1978 



Kucakova et al.|2011| . 

3fo7 fixed a, L K tends to infinity as the orbital radius approaches r = n co , 
risco (jjK(risco) = 2V3o^M and G r (ri sco ) = -^4- (l " 



The angular momentum L K as a function of r 
then monotonically decreases until it reaches its minimum value for r = 
where ri sco = 6M is the last stable circular orbit radius for a test particle in the Schwarzschild geometry. Finally it increases 
for r > ri S co- The angular momentum L K is a monotonically increasing function of a, and in the boundary a = it is L K = 0. 



3.1 Radial pressure gradient G r {L) vs angular momentum L 

Figure [T]-ng/i£ panel illustrates the sign of the radial pressure gradient G r as a function of the dimensionless angular momentum 
C = L/(Mcr) and distance from the attractor r. As we noticed in the previous Section, G r < (pressure decreasing) in the 
range 2M < r ^ ri co for every value of the angular momentum. When r > ri co , G r < (pressure decreasing) for < L ^ L K 
while G r > (pressure increasing) for L > L K . 

We are now interested in finding explicitly the critical points for the pressure, i.e. to find the solutions to Eq. [12] For a 
test particle within the effective potential V sc , the angular momentum L is a constant of motion. The particle motion is then 



2 The angular proper velocity of the fluid with L = is <E> K = L K /r 2 cr 2 . This function has no critical point. It is defined in r £]n co , 00 [, 
where it is a monotonically decreasing function of the orbital radius r/M, it increases approaching r = r^ co , and goes to zero at infinity. 

3 We restrict our analysis to a £ (0, 1] (the function is even in a £ [—1,1]) and L K ^ 



© 0000 RAS, MNRAS 000, [T]-?? 



On the Polish doughnut accretion disk 5 




Figure 1. (Color online) Left panel: r^ /M (black line) and r^ /M (orange line) as a function of L/M for different a. The dashed line 
marks = r^ = ri sco . This condition is fulfilled when C = 2\/3. Center panel: G r = constant as a function of t and r/M. The 
singularity at £ = £ r is marked with a thick black line. Dashed line marks r = ri co . Right panel: £k = L K /(Mcr) as a function of r/M. 
p' > (and p' < 0) in the white (gray) region. G r = when L = L K (solid line). Dashed line marks r = n sco , dotted line r = r\ co . 



<r=\ £(o-=l) 




Figure 2. Left panel: C on the equatorial plane a = 1 as a function of l/M G [0, 10] and r/M E]2, 10]. Central panel: C(a = 1) as a 
function of r/M, for different I. Asymptotically (r — > oo), it is L — > I. The dotted-dashed curve marks / = l r as a function of r/M: the 
angular momentum L is denned for I < l r . r\ co is a minimum of l r - Right panel: C{o = 1) = constant as a function of r/M and l/M. 



described by V' sc — (the prime ( ; ) stays for the derivative with respect to r) 
circular orbit radii for a test particle in Schwarzschild spacetime: 

(M ± VL 2 - 12cr 2 M 2 ) = Y ( M± gM ^ C? 12 ) 



which is equivalent to Eq. 



12 



2a 2 M 2 



Therefore, the 



(15) 



are solutions to Eq. 



12 



corresponds to the test particle stable orbit, and to the unstable orbit. The isobar surfaces 
(G r — 0) are therefore located at r = • Figure ^eft panel describes the behavior of rj^ as a function of L and a. We 
underline that the matter distribution around the accretor is not spherically symmetric, and hence it is not independent on 
a. A complete characterization of the equipotential surfaces as a function of a can be found in Appendix |A| . 



3.2 Angular momentum L vs. fluid angular momentum 



It is also possible to describe the motion of the fluid orbiting in the Schwarzschild background in terms of its angular momentum 
/, denned as: 

In fact, equipotential surfaces define the marginally stable configurations with respect to the axisymmetric perturbation 
( Seguin||1975 ), characterized by / constant. These configurations have been detailed studied in Kozlowski et al. (1978). From 
Eq. |16| and the definition of the effective potential in Eq. [7] the following relation holds: 



(17) 



e~ v r 2 G 2 l- 2 - 1 
The boundary case / = corresponds to L — 0. 

We detail the angular momentum L on the equatorial plane a = 1 as a function of / ^ and r in Figs. [2] The two 
definitions of angular momentum coincide in the asymptotic limit of flat spacetime (r — > oo). It is manifest from Fig. [2] 
central panel that it is L = I ^constant for sufficiently large distances from the accretor. L is not denned everywhere in the 
plane (r, I): from Eq. [TtI we notice that L exists for ^ / < / r , where: 



(r - 2M) 



(18) 
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is denned for r > 2M. In the boundary case r = 2M, L = 0. For fixed a, / r is increasing for r > n co , is minimum for r = n co , 
and is deceasing for 2M < r < ri co - It is also progressively larger when approaching the equatorial plane a = 1, where it is 
maximum. 

We now face the problem of finding a relation to link (r, a) and / to the condition of the existence of L, and therefore of 
the velocity of the fluid <I>. For this purpose we re-define the radii in Eq. [l5| 



r+/M 



2 VI 2 cos ^| arccos 



V3 



"/M : 



2^Fs 



Hi) 



V3 



(19) 



introducing the dimensionless quantity £ = l/(o~M) > (£ r = l r /(crM)). With this definition, rf 



n co for £ = 3^3. 



approaches the horizon r = 2M as the angular momentum £ increases (see Fig. [3] upper right panel). The angular momentum 
L and the velocity <E> are not defined inside the region [r~,r+] (see Figs. [3] bottom panel and upper left panel). The region 
[ r i "5 r t\ increases with £. It varies also for different a: its behavior as a function of a is detailed in Appendix [a| 

Now, we investigate the critical points of the angular momentum L as function of r, solutions of L' = 0. For r > n co , L 
is a constant with respect to the orbital radius when / = / K , where 



WM-„ M(r _ 2M)2 

is the Keplerian angular momentum of the fluid. The critical points of L are, thus: 



r-/M 
r + /M 



£ 2 -2^/£ 2 (£ 2 -12) i 



£ 2 + 2y/£ 2 [£ 2 - 12) cos 



7r + arccos 



2 (54 - \U 2 + ^ 4 
(£ 2 (£ 2 - 12)) 3/2 
2 (54 - 18^ 2 + £ 4 
[£ 2 (£ 2 - 12)] 3/2 



(20) 

(21) 
(22) 



where r_ ^ r+, and r_ = r+ = ri sco when £ = 3^3/2. 

The angular momentum L is a decreasing function of r/M in all r > n co with ^ K < ^ < £ r , it increases with r in 
2M < r < n co with < £ < £ r and in r > r Zco for < £ < £ K (see Figs. |3}. 



3.3 Radial pressure gradient G r (l) vs fluid angular momentum 

The radial gradient G r , as a function of £ = 1/ (jM), can be written as: 



G r — 



Mr 6 



2 M 2 (r-2M) 2 



(23) 



r(r - 2M) [r 3 - £ 2 M 2 (r - 2M)] 
for a / (see Figs. [TJ. It is not defined in £ = £ r = ^/r 3 /M 2 (r - 2M) and in r = rf. 

We studied the sign of G r as a function of L in Sect. |3.1| here we face the problem for G r as a function of £, the fluid 
constant of motion. In Appendix [A] we will discuss the G r sign in terms of (/, a, r) explicitly. According with the results found 
in Sect. 3.2 G r — (isobar surfaces) for the radii r± with Keplerian angular momentum £ K . Thus, G r < in: 

2M < r < nco for ^ £ < £ r and £ > £ K , 
r > rico for < £ < £ K and £ > £ r . 

In terms of the radii r± and rf, G r > in: 

0^£<3 v / 372 in r > 2M , 

^Sy^ in r>2M, r + r lsco , 

3 V / 372<^^3^3 in 2M < r < r_, r > r+ , 

£ > 3^3 in 2 < r < r z ~, r_ < r < 

These intervals are portrayed in Figs. [4] upper panels. 



r > r+ . 



In Sect. 3.2 we verified that L ; (/ K ) = 0, where L(l K ) = L K for r > ri co . Here we showed that Z K satisfies the condition 



Gt(Ik) = and therefore we can claim that l K is also a critical point for the pressure p. This is illustrated in Figs. [4] 



4 THE POLAR ANGULAR PRESSURE GRADIENT G# 

We now concentrate our attention on the polar angular pressure gradient From Eq. [To] 

© 0000 RAS, MNRAS 000, [l]-?? 



On the Polish doughnut accretion disk 7 






0</<3V3/2 


3V3/2 <f<3VJ 




2M<r<rf 


L' >0 


L'>0 




rf<r<r_ 


L' >0 


L'>0 




r_<r<rj + 


L' >0 


L'<0 




r, + <r<r + 


L' >0 


L' <0 




r>r + 


L' >0 


L'>0 





Figure 3. (Color online) Upper left panel r± and r l as a function of £. r+ = r_ = ri sco = 6M (dashed line) when £ = 3^3/2 
(dotted-dashed line). = = r\ co = 3M (dashed orange line) when £ = 3\/3 (dotted line). White region corresponds to L' > 0, light 
gray region to L' < 0. In the gray region the function L is not defined. Upper right panel 7y (black curve) and (orange curve) as a 
function of £ = l/(crM). r+ = r 



3M (dashed line) when £ = 3\/3 (dotted line). L is not denned in the gray region. Inset panel 
(black curve) and (orange curve) as a function of a for different values of I. Dashed line marks r = n co , where = . Lower panel 
table summarizing the intervals where L' > (light-gray) or L' < (gray). Black boxes correspond to the interval [r_,r+], where the 
function L is not denned. 



(24) 



P + V 

On the equatorial plane a = 1, it is = 0: the angular gradient of pressure on the equatorial plane is always zero. d$p = 
also for L = 0. This implies that the case L = leads to a zero polar gradient of p, while G r (L = 0) = —Mj ((r — 2M)r) < 0, 
i.e. the pressure decreases as 1/r 2 as approaching the horizon r = 2M. 

For L = L K , d r p = and therefore p is a function of & only. In general, for a / 1 it is \7#p(L — L K ) ^ 0. In fact it results: 

v^(Lk) _ mVT^ 



G$(L — Lk) 



3M)' 



(25) 



(see Figs. |5j|J The quantity in Eq. 25 is not denned in r — ri co and, as for L K , in the region r £ (2M,n co ]. This means 
that there is a fluid configuration with pressure constant along the orbital radius and a pressure variable from plane to plane 
extended at radial distance r > ri co . G$(L = L K ) admits critical points for a = 1. 



4.1 Pressure gradient ratio n r # vs angular momentum L 

We define the pressure gradient ratio II r ^ = n r #(L, r, a) as II r ^ = V r p/V^p. It is clearly II r ^ = for L = L K and it is not 
denned in a = 1. We rewrite this function as: 



n 



icottfn^Mi 



■3M)\ 



C 2 r(r - 2M) 
This ratio is equal to one for: 



(26) 



4 G$(L = Lk) is written as an odd symmetric function of cr, but we are in an even symmetric theory in cr, as we have not explicitly 
given any further constraints to the unknown functions (p,p), we implicitly take a > and y/l — a 2 ^ on? G (0, 7r/2]. 
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0</<3V3/2 


3V3/2 </<3V3 


/>3V3~ 


2M<r<rf 


p'<0 


p' <0 


p'<0 


rf <r<r_ 


p'<0 


p' <0 


p'>0 


r_<r<r 2 + 


p'<0 


p' >0 


p' <0 


rj*"<r<r + 


p' <0 


p' >0 


p'>0 


r>r+ 


p'<0 


p' <0 


p'<0 



Figure 4. (Color online) Upper left panel: r± and rf 1 as a function of £. White regions correspond to p' > and gray regions to p' < 0. 
Radii rf 1 are defined for £ > 3\f?>: the gradient G r is not defined in rf 1 . For 3^3/2 < £ < 3\f?> the radius r_ is a minimum of the 
pressure, and r+ is a maximum. For £ > 3^3 the radius r_ and r+ are both maximum of the pressure. r±, and are plotted as a 
function of £, r+ = r_ = n sco (dashed line) when ^ = 3y / 3/2 (dotted-dashed line). = r~[ = r\ co when £ = 3\/3 (dotted line). Inset: 
zoom of the region r/M £ [2,3]. Right panel: £ r (black line) and ^ K (orange line) as a function of r/M. £ r = £ K = 3\/3 in r = r\ co 
(dashed line), where G r = — (1/3)M _1 . The function is not defined in £ = £ r , while G r = in £ = £ K . In the gray regions p' < 0, in the 
white regions p' > . The angular momentum ^ K has a minimum in ri sco (dashed line). Lower panel: table summarizing the regions of 
increasing (jp' > 0-light-gray shaded regions) and decreasing (p f < 0-gray shaded regions) pressure p. 

^ G#(L=L K ) 
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1.0 



Figure 5. left panel: G$(L = L K ) as a function of r/M G]3, 20] for different values of a. The black dashed line marks r = ri sco . No 
solutions of dtfp = 0, and no test particle circular orbits exist for r/M £ (2, 3]. Right panel: Curves G$(L = L K ) ^constant as a function 
of r/M and a. 



r(r — M) — 3M 2 ' Tn ^ M 2 



C 2 



+ 



C? (13£ 2 - 12) 



C 2 -l ■ V (£2-1) 
The ranges where 11^ is larger, lower or equal to one are summarized in Fig. [H] left panel 



(27) 



4.2 Analysis of the pressure gradient ratio U r $ vs angular momentum / 

We consider now the ratio H r $ as function of /. We then define the function U l r $ — II^(r, £): 
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10 




Figure 6. (Color online) Left panel: as a function of C. Black line marks C = 1. In the light gray region 11^ < 1, in the white region 
11^ > 1. Right panel: r l n , as a function of t. Black line marks the value t = 1, dashed line t ~ 6.35, dotted-dashed line £ = 3y^3, 
dotted line r = | (l + -s/13) M. In the black region < r < the function II J,^ is not denned. In the gray region 11^ < 1, in the 



white region II 1 ^ > 1. In r = (black line) it is II 1 ^ = 1. 



n;^ = |cot^Mn r ^| 



|r 3 -fM(r - 2M) 2 



(28) 



^ 2 r(r - 2M) 2 

in the range of existence of L(l). Then we exclude the range r7 < r < r+. The solution of 11^ = 1 are r = r l n > M(l + >/l3) /2 



and £ = £ u G (M ( n-)), where £ (n -) = y/(41 + llv 7 !^) /2 = 6.35063 satisfies the condition r^(£ (n -)) = rj~(£( n _)) 
2.30278 M (see Fig.[6| riflffa pane/), with: 



(r + M)(r - 2M) 2 



and 



rk/Af: 



£ 2 



+ 2r 



(£ 2 



e 4 V (f 2 - 1) ; 



(1 - £ 2 ) (i> 4 - Al 2 - 2) 



(29) 



(30) 



The ranges where U l r # is larger, lower or equal to one are summarized in Fig. [H] right panel. 



5 THE BOYER POTENTIAL 



The disk fluid configuration in the Polish doughnut model has been widely studied by many authors (see for example Frank 
et al. 2002 Abramowicz& Fragile 2011). In particular, an analytic theory of equilibrium configurations of rotating perfect 

Boyer| Jl965| . The "Boyer's condition" states that the boundary of any stationary, 

constant. For a barotropic fluid the surfaces of 



fluid bodies was initially developed by 



barotropic, perfect fluid body has to be an equipotential surface W(l,$) 
constant pressure are given by the equipotential surfaces of the potential defined by the relation 

CPi n dp 



Jo 



-(W-W i7 



(31) 



/o P + P 

where the subscript "in" refers to the inner edge of the disc. It is important to notice here that, in the newtonian limit, the 
quantity W is equal to the total potential, i.e. to the sum of the gravitational and of the centrifugal effects. 

generic rotation law Q — are 



As mentioned in Sect. 3.2 all the main features of the equipotential surfaces for 
described by the equipotential surface of the simplest configuration with uniform distribution of the angular momentum 
density /, which are very important being marginally stable (Seg uin|1975 ). At the same time, the equipotential surfaces of the 



marginally stable configurations orbiting in a Schwarzschild spacetime are defined by the constant Z. It is therefore important 
to study the potential W = W(l, cr) and to compare it with W = W(L,a). 

We can classify the equipotential surfaces in three classes: closed, open, and with a cusp (self-crossing surfaces, which 
can be either closed or open). The closed equipotential surfaces determine stationary equilibrium configurations: the fluid can 
fill any closed surface. The open equipotential surfaces are important to model some dynamical situations, for example the 



formation of jets (see e.g. 


Kucakova et al.| 


2011 


|Stuchlfk et al.||2009| 


Lei et al.||2009; 


Abramowicz et al. 1998; 


Stuchlfk & 


Slany||2006| |Abramowicz 2009; Stuchlfk & ] 


<ovaf 2008; Rezzolla et al. 


2003; StuchKkl 2000; Slany & Stuchlfk 2005; Stuchlfk 


et al.||2000|). 
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1= constant L= constant 




t X 

Figure 7. (Color online) Left panel: (orange curve) and (black curve) as a function of £. The black region < r < r^~) is 
forbidden. V sc > 1 in the white region, V sc < 1 in the gray region. Dotted line marks £ = 4. Right panel: f^, (orange curve) and 
(black curve) as function of C. V sc > 1 in the white region, V sc < 1 in the gray region. Dotted line marks C = 4. 



The critical, self-crossing and closed equipotential surfaces are relevant in the theory of thick accretion disks, since the 
accretion onto the black hole can occur through the cusp of the equipotential surface. According to Paczyhski (lAbramowicz 



et al. 1978 Kozlowski et al.|[l978 Jaroszynski et al."|[l9 80i Abramowicz 1981), the accretion onto the source (black hole) is 
driven through the vicinity of the cusp due to a little overcoming of the critical equipotential surface W CU s P by the surface of 
the disk. The accretion is thus driven by a violation of the hydrostatic equilibrium, clearly ruling out the viscosity as a basis 
for accretion ( |Kozlowski et al!]|1978 ). In the Paczyhski mechanism the disk surface exceeds the critical equipotential surface 



Wcus P giving rise to a mechanical non-equilibrium process that allows the matter inflow into the black hole. In this accretion 
model the cusp of this equipotential surface corresponds to the inner edge of the disk. 
We calculate now the Boyer potential for our system integrating Eq. |11| 



/" 

J Pi: 



dp 
P + l 



dr 



The integration range [ri n ,r ou t] C [2M, oo] is the range of existence for V sc (r, /), (r 
condition implies that we are excluding the range [t^ ,r^~]. The general integral is: 

dp 



/; 



P + V 
where 

W = In 



W. 



(r - 2M)r 2 



r 3 -£ 2 M 2 (r - 2M) 



= hiV sc 



In particular it results: 
lim W — 0, lim W 

r— >oo i — >2M 



+ 00 



(32) 

2M) [r 3 - l 2 M 2 (r - 2M)] > 0. This 



(33) 



(34) 



(35) 



(see also Stuchlfk et al. 



2000 



Stuchlfk & Kovaf 



2008 



Abramowicz 



2005 



Abramowicz et al.|1980| |Rezzolla et aI.|2003) .~Clearly it is W = where V s 



Stuchlfk & Slany 



2006; 



Slany & Stuchlfk 



2005 



1, and W ^ where V sc ^ 1. Its maxima and 
minima are the same than the Schwarzschild effective potential, V sc = V sc (l,r). In particular, we are interested to study the 
equipotential surfaces, defined by the condition W = constant, that coincide with the surfaces V sc = K > 0, being K = e c the 
energy of a test particle circularly orbiting around the source. In fact the surfaces of constant Boyer potential determine the 
shape of the torus (disc). We study these surfaces as a function of L = constant and / = constant, respectively. 



5.1 Analysis of the Boyer potential vs the angular momentum L 



We consider the Boyer potential W(L,r) in Eq. 34 as function of the angular momentum L. The condition W = is satisfied 
in: 

f± = *t [C 2 ± ^C?{C? - 16)) (36) 
(see Fig. [7] right panel). 
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Figure 8. (Color online) ffc 1 , ffc 2 , ffc 3 as a function of £ 2 for different values of K in the ranges: [0,2^2/3] (left panel), [2^2/3,1] 



(central panel), K > 1 (right panel). Dotted lines are the momentum £ 2 and (C K ) 



±N2 



We notice that two relevant cases occur when the angular momentum is C 2 = C\. where Cf. 
< K < 1, W < 0, in 2M < r < f K , where f K ; (II) when K ^ 1, W ^ in r > 2M 



2 _ r 2 [2M+(K 2 -l)r] 
M 2 (r-2M) ' 



K4 



The solutions of the equation = ln(K) can be describer in terms of the energies: 

Ko 



54 

of the radii: 

f fcl /M = 

^ 3 /M = 
(see Figs, [s] 



1 / 36 ±A /^i^ + £ 2 



£ 2 



2 + 2zu sin [| (tt + 2 arccos -0)] 

3 (K 2 - 1) 
2 + 2zu cos [| (tt + arccos i/j)] 
3(K 2 - 1) ' 
— 2 + 2-ccj cos (| arccos^) 
3 (K 2 - 1) ' 

where: 
+ 9 (3K 4 - 5K 2 + 2) C 2 



1 

7! 



3£ 2 



C 2 



and of the angular momenta: 



w = (K 2 



1) 



4 + 3 (K 2 - 1) C? 



(K 2 - I) 2 



3(K 2 -1)' 1 Kj ~ 2 



27K 



4 36K 2 + j 



K 2 



K 2 (9K 2 - 8) 3 



(K 2 



These solutions are summarized in Figs. [9] 

The critical points of the angular momentum C 2 are: 



rt/M : 



K 2 (9K 2 



3K 2 



(K 2 



K 2 - 1 



(I) when 

(37) 

(38) 
(39) 
(40) 

(41) 
(42) 



(43) 



When 2y / 2/3 < K < 1 there are two critical points fj, while when K > 1 there is only r = fjj" (see Fig. [9] upper right panel). 
For K = 1 it is fjf = r m b , where r m b = 4M is the marginally bounded orbit for a test particle in the Schwarzschild spacetime, 
and in K = 2y / 2/3 it is fjf — ri sco . 

Closed surfaces: the conditions for the existence of closed surfaces can be obtained by noting that, in cartesian coordinate 
(x,y), the closed surfaces should satisfy the condition V sc (x = 0) = K with three solutions, say y = {2/1,2/2,2/3}- The closed 
surfaces then exist when 2^/2/3 < K < 1 and (C^) 2 < C 2 < (C^) 2 (see Fig. [9] lower panel), with (C^) 2 ^ 16. In cartesian 
coordinate (x, y) the surfaces V sc = K are: 



(2L 2 + 2y 2 



(L 2 + y 2 -y 2 K 2 y 

The maximum diameter (x = 0) of the closed Boyer surface lies between the points y = 2/2 and y = 2/3 , where 



2/i /M 



q — 60 sin 



[tt + 2 arccos e) 



(44) 



(45) 
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Figure 9. (Color online) Upper left panel: K+ (black curve), K~ (orange curve) and K a (dashed thick curve) as a function of C? . The 
regions of existence for the radii , r^ 2 , rfc 3 , solution of W = ln(K) are marked in white, gray, light-gray. Black regions are forbidden. 
Dotted-dashed line marks K = 2 a/2/3. Upper right panel: fx as a function of K. Gray curve is fi", dashed curve r^. In the white region 
C 2 = C 2 ^ solution of W = ln(K). No solution exists in the black region. Lower panel: table summarizing the regions of existence of the 
radii r^ 1 , ffc 2 , ffc 3 . Black boxes are forbidden. 



2/2 /M 

2/3 /M 
being, 



^ — 60 cos 



- (tt + arccos e) 



q + 60 cos ( - arccos e 



16 + 8 (3K 4 - 4K 2 + 1) (L/M) 2 + (K 2 - l) 2 (L/M) 4 



(K 2 - l) 4 

64 + 48 (3K 4 - 4K 2 + l) (L/M) 2 + 6 (K 2 - l) 2 (9K 4 - 6K 2 + 2) (L/M) 4 - (K 2 - l) 3 (L/M) 6 



(k 2 - iyo 3 



6 [2+ (K 2 - 1) (L/M) 2 



(46) 
(47) 

(48) 
(49) 
(50) 



(k 2 - iy 

(see Fig. [To] left panel). 

Fig. |10| right panel portrays a closed Boyer surface. We can characterize this surface introducing the following parameters 

(i) the surface maximum diameter: A = 2/3 — 2/2, 

(ii) the distance from the source denned as: S = y<2 — 2M, 

(iii) the distance from the inner surface: 5 = 1/2 — yi, 

(iv) the distance of the inner surface from the horizon: 8 = y\ — 2M, 

(v) the surface maximum height denned as: h = 2xm, 

(vi) the quantity Xu = yM- 
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Figure 10. (Color online) Left panel: y\ (black curve), 1/2 (orange curve) and 7/3 (dashed curve) as a function of L 2 for K = 0.96. Right 
panel: the closed Boyer surface at K = 0.98 and L 2 = 16M 2 . A = 1/3 — 1/2 is the surface maximum diameter, 6 = 1/2— 2M is the distance 
from the source, 5 = 1/2 — yi is the distance from the inner surface, 5 = yi — 2M is the distance of the inner surface from the horizon, 
h = 2x M and X h = jjm- 



(xm,Vm) is the critical point of the surface in Eq. 44 In what follows we find the constraints for the set of parameters 
{A, 5j 5) <5, hj Xh}. The point y% varies in the range 2M < y% < y s , and L 2 lies on the surface L 2 , which is a solution of V s = K 
on the plane x = for y > AM and K ^ 1, and for 4M < y < y s for K < 1: 



Vs 



2M 



IK 2 



2M 



K 4 y 6 



( y 2 _4M 2 )^ 



+ 1T 



K 2 y 2 



y 



4M 2 



(51) 



(see Figs. 11 ). A increases with the energy K, but decreases with the fluid angular momentum L. On the contrary, the distance 
from the source 5 = 2/2 — 2M increases with L and decreases with K (see Figs. 11 ). 
The maximum vertical distance of the closed surfaces is: 



y M /M 



3(L/M) 2 +4V6(K 2 -1) 



K 2 (L/M) 2 
" (K2-l) 3 



(K 2 



K 2 (L/M) 2 
" (K 2 -l) 3 



3(K 2 



(52) 



As shown in Figs. 



11 



the maximum xm — x(^/m), and consequently the height h, increases with K and decreases with L 
On the contrary yM increases with K and with L 2 , until it reaches a maximum and then decreases. The distance 5 = \ 
increases with L 2 and decreases with K. The distance 8 = yi — 2M increases with the energy K and decreases with L 



yi 



Cusps: the cusps, i.e. the self-crossing surfaces W = constant, correspond to the maxima of the effective potential as a 
function of r: open surfaces are maxima with energy K > 1, and closed surfaces are maxima with energy K < 1, as outlined 



in Fig. 12 upper. It is therefore important to consider the function V' in the regions of closed and open Boyer surfaces. From 
Sect. 3.1 we know that the solutions of the equations V' — that correspond to maxima of the effective potential are located 



on the radius r = r L , in particular there are closed surfaces with a cusps in r\ 



located in r T 



cusp (see Fig. 



12 



upper). 



and 12 ^ C 2 < 16 , viceversa maxima 



- l > r mbo with C ^ 16 are open surfaces with 

Remark: integrating Eq. [32] with V sc as function of the constant of motion L, we obtain the following expression for the 
potential: 



W : 



2£ 4 M(r + 3M) + (4 + £ 4 ) r 2 [ln(r/M - 2) - ln(r/M)] 



8r 2 



(53) 



The solutions of W = c ensure the existence of closed, open and self-crossing surfaces. 



5.2 Analysis of the Boyer potential vs the angular momentum / 

We face now the analysis of the Boyer potential in terms of the fluid angular momentum / = constant. Firstly, we note that 
the Boyer surface in Eq. 



34 



is not denned in the region [r z - ,ry ]. We detail the study of the sign of V sc (l,r) in Figs. 13 (see 
also Fig. [7] left panel). With respect to the case of L — constant (see Figs.|9]and Eqs. 37 42), here new definitions have been 
introduced for the energies: 



K a 



' £ 2 — 36 — 2/3i sin arcsinai) 



3^ 2 - 81 



(54) 
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Figure 11. Upper left panel: y s as a function of the energy 2 a/2/3 < K < 1. Inset: distance of the surface from the r = 2M S = 1/2 — 2M 
as a function of L 2 , for different values of K. It increases with L 2 and decreases with K. Upper right panel: surface maximum diameter 
X = ys — z/2 as a function of L 2 for different values of K. Lower left panel: maximum xm = x(um) as a function of L 2 for different 
energies K. Inset: yM as a function of L 2 , for different values of K. Lower right panel: distance 5 = y\ — 2M as a function of L 2 , for 
different values of K. Inset: distance 5 = y2 — yi as a function of L 2 , for different values of K. 



K h 



K 



with 



I 2 — 36 + 2/?i cos ( \ arccosai) 



U 2 - 81 



1 / Wl 2 (3l 2 -16) 



(2 3 3 9 ) - (108^2) ^ 2 + (6V51) 2 ^ 22 - 72* 23 + 



^3 , ^24 



— ^ 2 



£ 2 £? 



angular momenta: 
4 



a ~ 3K 2 (K 2 - 1 
and radii: 



-. (£) 



±\2 _ 



(9K 2 



+ 



(r - 27) 



27K 4 - 36K 2 + ; 



72 + (£ 2 - 24) 2 
(^ 2 - 27) 2 ' 



K 2 (K 2 -1) 2 K 2 (K 2 -1) 



r kl /M 

rk 2 /M 

rkjM 
where: 



2 [l + vjk cos (| arccos^fe)] 
3(K 2 - 1) ' 
2 (-1 + wk cos [§ arccos (-^fc)]) 



3(K 2 - 1) 



2 [-1 + 



zuk sin I o arcsm 



in^fc)] 



3(K 2 



(55) 
(56) 

(57) 

(58) 

(59) 
(60) 
(61) 
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Figure 12. Left panels: Effective potential V SC (L, r) as a function of r/M, for selected values of C 2 = L 2 /M 2 a 2 . Dashed lines mark rj . 
Inset: the self-crossing Boyes surfaces corresponding to the maxima of the effective potential as function of r: maximum with energy 
K > 1 are open surfaces and maximum with energy K < 1 are closed surfaces. Right panels: Plots of the effective potential V sc (l,r) as 
function of r/M, for selected values of t 2 = l 2 /M 2 a 2 . Dashed lines are r±. Correspondingly the self-crossing Boyer surfaces are plotted 
in the inner panel. They correspond to the maxima of the effective potential as function of r: maximum with energy K > 1 are open 
surfaces and maximum with energy K < 1 are closed surfaces. 



8 + 9K 2 (3K 4 - 5K 2 + 2) f 2 / 4 + 3K^ (K^ - 1)P 
*u = ^ , WfcS (K-l)^ . (62) 

The surfaces W(l,r) — constant exist in all the spacetime with angular momentum ^ 2 = £l and energy K = Kk, where: 



ik " K 2 M 2 (r-2M/ ' < M 2 (r-2M)' Kk = T \J r 3 -W(r-2M)' ^ 
In particular in the case < K < 1 it is 2M < r ^ 2M/(1 - K 2 ). 
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Figure 13. (Color online) Upper panel: K a (black curve), (orange curve), KJ (dashed thick black curve) and (dashed thick 
orange curve) as a function of £ 2 . The colored regions mark the existence for the radii , r^ 2 , rj~ 3 • The black region is forbidden. Dotted 
lines are £ 2 = 16 and £ 2 = 27/2. Lower panel: table summarizing the regions of existence for the radii r^ 1 , r^ 2 , r^ 3 ■ Black boxes are 
forbidden. 



Closed surfaces: following the same procedure outlined in the previous Subsection for V SC (L, r), we find that closed surfaces 
of the Boyer potential in the cartesian coordinate (x,y) are in the regions 2\f2/3 < K < 1 and ) 2 < £ 2 < (^) 2 , where 
£ 2 k > 27/2 (see Figs." - 



13). The surfaces are therefore: 



x = ± 



2M (K 2 1 2 + y 2 



K 2 (l 2 - y 2 ) + y 2 



(64) 



(see Fig. 14 panel). 

The maximum diameter (x = 0), is denned by the points y = y2 and y = 2/3, where 



yi/M 

y 2 /M 

ys/M 



S k - 60 k sin 



(7r + 2 arccos£fc) 



^ - 6^fc cos 



-(tt + arccos^fc) 



^fc + 60 k cos 



■ arccos 



(65) 
(66) 
(67) 
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Figure 14. (Color online) Left panel: y±/M (black curve), yi/M (orange curve) and yz/M (dashed curve) as a function of l 2 /M 2 for 
K = 0.96. Right panel: the closed Boyer surface at K = 0.98 and I 2 = 16M 2 . A = ys — y2 is the surface maximum diameter, 6 = y2 — 2M 
is the distance from the source, 5 = y2 — yi is the distance from the inner surface, 5 = y± — 2M is the distance of the inner surface from 
the horizon, h = 2xm and Xh = yyi- 



see also Eqs. 45 47| , being 



16 + 8K 2 (1/M) 2 + K 4 [(l/M) 2 - 32](1/M) 2 - 2K 6 [(1/M) 2 - 12](1/M) 2 + K 8 (l/M) 4 



(K 2 - l) 4 



64 + K 2 (l/M) 2 (48 + K 2 [48 (3K 2 - 4) + 6 (K 2 - l) 2 (9K 4 - 6K 2 + 2) (l/M) 2 - K 2 {K 2 - l) 3 (//M) 4 j ^ 



(K 2 -l) 6 ^ 



6[K 4 (l/M) 2 -K 2 (l/M) 2 + 2] 
(K 2 - l) 2 



(68) 

(69) 
(70) 



(see Fig. 14 left panel). 

The maximum height for the surface is: 



Vm/M 



\ 



3K 2 (1/M) 2 + 4^6(K 2 - 1) 



K 4 (l/M) 2 



(K 2 



(K 2 



K 4 (l/M) 2 
" (K2-l) 3 



(71) 



3(K 2 -1) 

The surfaces for / = constant are larger then those for L — constant. The two cases are compared in Figs. [l5| the maximum 
diameter \(l) > A(L); the distance from the source 8(1) > S(L); the distance of the inner surface from the horizon 5(1) > S(L); 
its maximum height h(l) > h(L) and finally the quantity Xh(l) > xh(L)\ the distance from the inner surface 8(1) < S(L). 

Cusps: the closed (open) self-crossing surfaces W — constant are located on the maxima of the effective potential with 



energy K < 1 (K ^ 1) as function of r (see Figs. 12 bottom). The critical points are located on r±: the maximum is r_ when 
3^3/2 < £ < 3\/3. The energy V sc (r-) at the maximum r_ and at t — 4 is V sc (r-) — 1 and for £ < 4 it is V sc (t-) < 1. 



6 THE POLYTROPIC EQUATION OF STATE 

We consider the particular case of a polytropic equation of state p(r) — /cp(r) 7 , where the constant 7 is the polytropic index 
and k > is a constant. Using this relation in Eq. |3l]we have: 

^C t w^) = - {WM - Winn)) - (72) 

Integrating Eq. [72] we obtain: 



In 



(p (fc + p 1 " 7 )" 1 ^ 



= -W, for 7^1, (73) 
and 

In [((l + /c)p)^r] = -W, for 7 = 1, (74) 

(isothermal case). 

Solving Eq. [73] and [74] for p and using Eq. [34] we find respectively: 
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Figure 15. (Color online) Upper left panel: comparison of the closed Boyer surfaces with energy K = 0.98 and angular momentum 
I 2 = 16M 2 (black curves) and L 2 = 16M 2 (orange curve). It results x(l 2 = 16M 2 ) > x(L 2 = 16M 2 ). Upper right panel: comparison 
of the surface maximum point i/m/M. Inset: comparison of the surface maximum point Lower left panel: comparison of the 

distance of the inner surface from the horizon 6 = y± — 2M. Inset: comparison of the surface maximum diameter A = 2/3 — 2/2- Lower 
right panel: comparison of the surface distance from the horizon 6 = 1/2 — 2M. Inset: comparison of the distance from the inner surface 

S = V2 - yi ■ 



1 

(h = [k 
and 

Pk = V sc 



V sc 



(-1+7) 



for 7/I 



1 + 



for 7=1. 



(75) 



(76) 



In the following we adopt the normalization: p 1 = k 1 ^ 1 ^ p 7 , which is independent from k. The following limits are satisfied: 

lim p,y — 0, lim 0^ — 00, (77) 
and 

lim pk = 1, lim pk — 00. (78) 

r— >oo r — >2M 

We underline that for 7 = 1 it is p° k ut = pT (V s ° c ut /V^)'^ , and for 7 ^ 1 it is k{p)Z^ = (Hp)]' 1 + 1) (V™* /V™) ^ - 
1. In the polytropic case it is p' — p' /(Aryp 7-1 ), thus it is p' = when p' = and, being 7 > 0, the maxima (minima) of p 
correspond to maxima (minima) of p. 



6.1 The case: 7 / 1 

If the polytropic index is 7 7^ 1, the density p — p 1 is: 

with C = (V~ 2 )^ - 1 and (V- 2 ) = (^ - ^) > 0. 



(79) 
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Figure 16. (Color online) Left panel: j q as a function of q, where \q\ ^ 1 are integers. For large values of \q\, the index 7 tends 
asymptotically to 7 = 1. Gray region is forbidden. The maximum value of j q is j q = 3/2, the minimum is 79 = 1/2. C < for V S 2 C < 1 
with < 7q < 1 in Rj and V^ c > 1 with 7^ > 1 in Ru. Right panel: (black line), rj~~ (orange line), (thick black line) and 
(thick orange line) as a function of t. Dashed line marks t = 4, dotted line t = 3a/3, dotted-dashed line r = r m b , thick-dashed line 
r — r lco- White region correspond to the range Rj, gray region to Ru. 



We distinguish between two cases: 

(i) C > and the density p 1 is denned for all 7; 

(ii) C < and the density p 1 is denned for 7 = 7^ = 1 + where \q\ ^ 1 are integers (see Fig. 

The condition C > is satisfied in two cases: where Vg C < 1 and 7 > 1, in the ranges: 

^_fo<^<4 in r >2M 
1 _ I £ ^ 4 in 2M < r < r w r > r+ 

and where Vg C > 1 and < 7 < 1, in the ranges: 



16 



left panel). 



4<£< 3V3 in 
R u = { £ = 3V3 in 
I > 3^3 in 



< r < Ty 

< r < ri 



< r < r, 



r / Tic 



with 



16, 



(80) 



(81) 



(82) 



(see Fig. 16/ right panel). 

We can summarize as follows: when the polytropic index 7 = j q , the fluid density p is denned for the conditions Ri Ufti, 
when 7 / 7^ it is denned only for the conditions R T . In the following subsections we will discuss an example of 7 = 5/3 / 7 g 
and the particular case j q (q = 1) = 3/2. 



T/ie adiabatic case: 7 = 5/3 

We consider now the particular case 7 = 5/3. This polytropic index is adopted to describe a large variety of matter models, 
as the generic degenerate matter like star cores of white dwarfs (see e.g. Horedt 2004). 
The density p 1 is then: 



P5/3 



M 



2 2 \ 1/5 



2M 



1 



3/2 



denned in the range Ri as in Eq. [80] 

<£ <4 in r > 2M, 

t ^ 4 in 2M<r<v, r>4 

where lim Mr ± p 5/3 = 0, (see Figs. 



17) 



The critical points of p 5 / 3 can be found as solutions of p' 5 , 3 = 0: 
■■3y/3/2 in r = risco, 3y^3pl < i < 4 in r = r±, 4 in r = r+, 



(83) 



(84) 
(85) 



(86) 
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Figure 17. (Color online) Left panel: p 5 / 3 as a function of r/M for different values of £. Right panel: p 5 / 3 as a function of r/M and £. 
P5/3 is not defined in [r^,rjjr]. 



r=5/3 




p 5 / 3 (l/M=8) 



3.0 3.5 4. 





















■'■•Si 






- y \ 






-20 -10 




rsini?/M 


10 20 




Figure 18. (Color online) Left panel: and r± as a function of £. In the gray (white) region p' 5 / 3 < (p' 5 / 3 > 0), in the black region 
p' 5 / 3 is not defined. Central panel: p 5 / 3 = constant in the plane (x,y). Right panel: p 5 / 3 as a function of (x,y) for I = 8M. The black 
surface is \J x 1 + y 2 = 2M. 



and it is p' 5 / 3 > 0, (density increasing with the orbital radius) for: 
3v / 3/2<^^4 in r_ < r < r+, ^>4 in 7^ < r < r+. 



(87) 



We thus conclude that r_ is a minimum and r+ is a maximum of p 5 / 3 (see Figs. 18). 
6.1.2 The case: 7 = 3/2 

We consider the particular case j q (q = 1) = 3/2, which is in the extreme cases 1 < 7 < 5/3. 
The density p 3 / 2 / is then: 

1 2 



P3/2 



A I 



2 2 \ 1/6 



1 



K r-2M r 2 
defined in the range Ri U R u as in Eqs. |80pl 

<£<4 in r> 2M, 
4 < ^ < 3^3 in r > 2M, r^^ 

O 3^3 in 2M<r<r l ", r > 
where lim Mr ± p 3/2 = 0, (see Figs. 



(89) 
(90) 
(91) 



19) 



The critical points of p 3 / 2 can be found as solutions of p 3 / 2 = 0, or noting that p = \p> ' y/p. We summarize these results 
concluding that r_ is a minimum and r+ is a maximum of ps/ 2 - 



6.2 The isothermal case: 7=1 



For an isothermal equation of state (7 = 1), the solution of Eq. 74 is p = pk (see Eq. 75). This function is not defined in the 
range r~[~ ^ r ^ rf and the fluid angular momentum is < £ < £ r . 
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In order to describe the regions of maximum and minimum density, we study the function p' k . Being p = kp, it is clearly 
p — kp . Thus it is p — when p — 0. Moreover, being k > 0, the maxima (minima) of p correspond to maxima (minima) of 
p. The existence of critical points for the isothermal case is therefore studied in Sect. |3.3]in terms of the critical points of G r : 



£ = 3^/3/2 in r = n sco , 3^/3/2 < £ < 3^3 in r = r±, £ ^ 3V3 in r = r+, (92) 
or also r > ri co i — i^- 

7 THE FLUID PROPER ANGULAR VELOCITY 

This Section concerns with the analysis of the fluid velocity field, in particular we are interested in assessing the orbits and 
the plans where the fluid proper velocity is maximum or minimum. 
The fluid four velocity along the ip angular direction is: 



* = — = —J ° 2l2 (r-2M)r> (93) 

r 2 a 2 a 2 r 2 y a 2 r 3 _/2( r _2M) ' ^ ' 

where we are always considering L > and $ > We redefine $ as a dimensionless quantity: 

$a 2 _ M 2 / £ 2 (r-2M)r 2 ^ g4 ^ 



M r 2 y r 3 -£ 2 M 2 (r-2M) ' 

Clearly it is $ ^ when L ^ 0, and they have the same existence conditions (0 < / < l r , 2M < r < and r > r^). The 
behavior of the angular velocity as a function of a, /, and r are portrayed in Figs. |20| 

In order to derive the critical points of the proper angular velocity as a function of r, we consider the solutions of the 
equation & = 0. Apart from the trivial case / = and $ = 0, the proper angular velocity is constant with respect to the 
radial coordinate when: 



l = l ^flSi^ fOT 5 /2M<r<3M, (95) 
where < l x < 3V3M. In terms of the angular momentum and the radius, the orbits of $ belong to the planes: 



(r - 2M) 2 

for r x < r < ri co (cr — 1 when r = r x ) , (96) 



r 3 (2r - 5M) 



where: 



Tx/M ~ 8 + 2V8 + 3 6a 6 T T 8^3' 1 ' 



with 

1 1/3 



a = [-£ 2 [t - 54£ 2 + 1350 - 6^48^ - 2754£ 2 + 50625^] , (98) 

P = Jl6£ 2 + 75+- 2 ^ 2 ~ 36 ) + 8 a . (99) 
V a 

5 The convention adopted is that the matter rotates in the positive direction of the azimuthal coordinate cp. 
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Figure 20. Upper left panel: <J> = aU^/M on the equatorial plane a = 1 as a function of r/M and l/M. The function is defined for 
I G [0, l r ]. Upper right panel: <I> = constant as function of l/M and r/M on the equatorial plane a = 1. The curves have minima in 
r = r x . Lower left panel: <£> = constant in the plane cr and r/M for Z/M = 1. The curves have maxima in r = 5/2M (dashed line). Lower 
right panel: <E> = constant in the plane a and t for r/M = 5/2. No motion is allowed for t > i r (I = £ r , dashed line), and it is <E> ; = on 



As shown in Figs. 21 2.5M < r x < 3M and the orbital radius r x increases with £ — I /(Ma). An alternative analysis of the 
proper fluid angular velocity, as function of (Z,r, cr), is presented in Appendix [5] 



8 COMPARING THE FLUID RELATIVISTIC ANGULAR VELOCITY AND THE KEPLERIAN 
ANGULAR MOMENTUM 

In this Section we compare the fluid configurations with d r p / with the model p ^constant where the disk is geodetic by 
studying the relativistic angular velocity and comparing our results with the Keplerian definitions. In Sect. |3.2| we already 
explored the fluid behavior with respect to Ik- We summarize here the results considering the ratio: 



n LK L/L K $/$ K 



PM(r -3M)(r -2M) 



£ 2 M 2 (r-2M) 



(100) 



which is defined for < £ ^ 3\^3 in r > ri co , and for £ > 3\^3 in r > r+. The regions where L(l) is larger or smaller 



\im r ^ rico n L K = 0. In general, L > L K for fluid angular momentum / sufficiently large, for / sufficiently low or 



than L K = L(l K ) for different values of the fluid constant of motion / are portrayed in Figs. 22 upper left panel. We have 
linv^oo ri| K 

far from the source it is L < L K (see Figs 

We are now interested to define the critical points of the ratio IIlk- The solutions of the equation IIlk — are, in terms 
of the angular momentum: 



22 



upper left panel). 



_ r v / (10M -r)r - 18M 2 
Ulk = M(r - 2M) ' 

and in terms of the orbital radius: 



(101) 
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Figure 21. (Color online)Le/£ panel: r^ ^constant as a function of < a < 1 and < l/M < 3\/3. Inset r^ (black line) and (orange 
line) as increasing functions of £ = l/{aM). Right panel: r^/M (orange line) and r^/M (black line) as a function of a. Gray region 
corresponds to r G [ri C o, °°]- Dotted lines mark a = o\ for the selected l/M. 
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Figure 22. Upper left panel: r± and r^~ as a function of £. IIlk is not defined in the black region. IIlk > 1 in the white region and 
n LK < 1 in the gray region. Dashed line marks t = 3^3/2, dotted line t = 3a/3- Upper right panel: ru LK , rf, r± as a function of 
£. IIlk is not defined in r < (black region). For r < ru LK it is IIlk > (white region). The gray colored region corresponds to 
H LK < 0- When r = r± it is il LK = 1. Dashed thin line marks £ = 3 a/3- Lower left panel: table summarizing the regions IIlk > 1 
(light-gray regions) and IIlk < 1 (gray regions) for different values of r/M and £. IIlk is not defined in the black boxes. Lower right 
panel : n LK (^n iK ) as a function of r/M £ [3.76,7.88]. Dashed line marks IIlk = 1. Inset: n LK (rn LK ) as a function of £ E [0, 3-\/3]- 
Dashed line marks IIlk = 1. 
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rn LK /M = 



15 + V3\ 78 -U 2 



(£ 2 - 18) 2 



6^3(^-35) 



S-.V3 _ " v " ^ "± + ^ , (102) 

(103) 
(104) 

These are defined for n co < tu lk < (5 + y/7) and < lu LK < 3V3. The sign of IIlk as a function of r and I is summarized 
in Fig. |22[ upper right panel. In particular, it is manifest that rn LK is a maximum of n LK . The values of this maximum, 
Hu LK (£) = n LK (^ = £n LK ) and Uu LK (r) = n LK (r = rn LK ), are portrayed in Fig. 22 lower right panel. We note in particular 
that Hlk is an increasing function of I: this implies that the larger the angular momentum of fluid, the larger is the ratio 
between L and the angular momentum of the free particle Lk> 



with 



[54 (108 + 10^ 2 - t) + t + 12y/6y/-£ 2 (£ 2 - 27) (£ 4 - 54£ 2 + 756) 



1/3 



£ 4 - 2£ 2 (1S + k) + (12 + k)(27 + k) 



8.1 The fluid relativistic angular velocity £l: the Von Zeipel surfaces 



The fluid relativistic angular velocity £l is: 

n = £ = = ; ( r - 2M ) = i 

S 9<*p<~p r^o~ 2 I 2 

The function Q = cr^(^): 
_ ^M(r - 2M) 



is defined in r > 2M, and has a maximum in r = n co (see Fig. 23 Ze/fc panel). 

The surfaces known as the von Zeipel 's cylinders, are defined by the conditions: 



(105) 



(106) 



constant and Q = constant (see 



for example Abramowicz 1971 Chakrabarti 1991 1990). In the static spacetimes the family of von Zeipel's surfaces does 



not depend on the particular rotation law of the fluid, Q = Q(Z), in the sense that it does not depend on nothing but the 
background spacetime. In the case of a barotropic fluid, the von Zeipel's theorem guarantees that the surfaces Q = constant 
coincide with the surface with constant angular momentum. More precisely, the von Zeipel condition states: the surfaces 
at constant pressure coincide with the surfaces of constant density (i.e. the isobar surfaces are also isocore) if and only if 
the surfaces with the angular momentum I = constant coincide with the surfaces with constant angular relativistic velocity 
( |Kozlowski et al.|1978| |Jaroszynski et al.|1980| |Abramowicz]|1971| |Chakrabarti||1991[|1990 ). 
The surfaces Q/l = s > 0, being s a constant, are defined by: 



~/M 



7T + arccos 



3V3 



V3 



rf/M 



3V3 



V3 



(107) 



These are a function of the product cr 2 s, and they exist for < o~ 2 s ^ (1/27)M 2 . In particular, a 2 s = (1/27)M 2 corresponds 



to r s 



: Tic 



see Fig. 24 upper left panel). In the (x, y) coordinates, the £1/1 = s > surfaces read: 



2M 5 



M 4 ■ 



sy z 



(see Figs. 24). 



We are now interested in characterizing the angular velocity with respect to the Keplerian velocity Q K 
therefore consider the dimensionless difference An = (ft — f2 K )Mcr: 



r 3/2 



Mcr 



^M(r - 2M) - Vr' 3 M 



(108) 
We 
(109) 

U L - 

An = for r 00, A^(r = 2M) = -1/2^2, A (r Zco ) = (£ - 3^3) /27, and A n (r Zsco ) = (2£ - 3^6) /108 (see Fig.|3 
right panel). An(r — 2M) is negative irrespective of the value of £. The sign of An(ri co ) and An(n sco ), on the contrary, 
depends on the fluid angular momentum. The sign of Aq as a function of r and £ is summarized in Figs. [25] 

The difference Aq is maximum in r > ri co when £ — (3/4)^ K - Notably, £ K is not a critical point of the angular velocity 
difference An. An always increases in 2M < r ^ ri co for all £, and for < £ < (3/4)^ K in the region r > n co . Other critical 
points are for £ ^ 27/8 in r = r^~: 

To/M = ^(,2£ 2 -V2^PW 



81) cos 



(tt + arccos 6) 



(110) 
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Figure 23. (Color online) Upper left panel: the relativistic angular velocity Q in units of 1/M, for cr = 1 as a function of r/M, for 
different values of £. Dashed line marks the Keplerian velocity f2 K - Dotted line marks r = ri co = 3M, a maximum of Q. For £ = 3^3/2, 
it is Q = Qk in r = r\ co = 3M; for £ = 3^/3, it is Q = Qk in r = ri sco = 6M and in r = 6(2 + y/3)M. Upper right panel: Aq as a 
function of r/M, for different values of £. Dashed line marks the maximum of Aq. Lower left panel: rj and r± as a function of £. In the 



white (gray) regions it is A^ > (A^ < 0). Aq = when r = r±, and Aq is maximum in r : 
panel: Aq^) as a function of I. Inset: A^frj") as a function of £. 



and minimum in r : 



Lower right 



rt/M 
with 



27 



2£ 2 + V2y/£ 2 (8£ 2 - 81) cos 



■ 7776r + 1 



16>/2 (-8£ 2 -81)] 



3/2 



(ill) 



(112) 



where r 



9M for 



r > r„ , while it decreases in 



27/8. Aq always increases in < £ < 27/8, and for £ ^ 27/8 in the regions (2M,r~) and 
^ 27/8 in the region (r~,r^"). The radius is a minimum and r~ a maximum point of An. 
are portrayed in Figs. 



The critical points values An (7*0 ) are portrayed in Figs. 23 lower panels. 

We analyze the profile of the proper angular velocity, <£, and the relativistic angular velocity, respect to the corre- 
sponding Keplerian quantities for the case of a isothermal matter. Clearly, with respect to the general analysis, here we should 
take into account the specific density and pressure profile emerging by the choice of the particular equation of state, p = kp 1 . 



9 SUMMARY AND CONCLUSIONS 



The analysis of a stationary axisymmetric configuration of material, in equilibrium in a Schwarzschild spacetime, as emerging in 
the Polish doughnut framework constitutes a timely question in view of the interest in astrophysical sources, possibly resulting 
from super-Eddington accretion onto very compact objects, like Gamma-Ray Bursts, Active Galactic Nuclei, binary systems 
and Ultraluminous X-ray sources (see e.g. Fender & Belloni 2004; Soria 2007). The most characterizing features of the Polish 
doughnut approach is the thickness of the matter distribution across the equilibrium and the existence of a region enveloping 
the horizon surface, where the fluid can, in principle, infall onto the black hole. The first property is of impact for a comparison 
with the wide spectrum of numerical simulation of a thick disk (see Straub et al.|2012 Font 1 2 003 Abramowiczfe Fragile|20TT 
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Figure 24. (Color online) Upper left panel: rJ/M as a function of (cr 2 s). Dashed lines mark the conditions r = ri ccn cr 2 s = (1/27)M 2 . 
The white region (r < ri co ) corresponds to Q'/l > 0, the light-gray region (r > n co ) to Q' /I < 0. Q'/l = in r = r\ co . Upper right panel'. 
surfaces Q/l = s with r/M = constant in coordinate x = rcos#, y = rsin#. The surface s = (1/27)M 2 is marked with a number. The 
gray region marks the circle r ^ 2M. Arrows follow the increasing (decreasing) direction respect of s. Lower panels: surfaces Q/l = s in 
units of M 2 as a function of x = rcos$, y = rsini? in units of M, for different values of s. The gray region underlines the circle r ^ 2M. 
Dashed circle corresponds to r = ri co . 



for recent examples). The non- negligible depth of the accreting profile is typical of the regime where the gravitational effects 
are strong and it takes a relevant role in all those extreme phenomena associated with the gravitational collapse, characterized 
by a violent energy-matter release from the central compact object. The second aspect is relevant because the Polish doughnut 
model can account for a non-zero accretion rate of the torus even when the dissipative effects are negligible. This is in contrast 
to the original idea by Shakura &; Sunyaev ( 1973 ) that the angular momentum transport is always allowed by the shear 
viscosity of the accreting material. Indeed, the accreting plasma is in general quasi-ideal and the emergence of dissipative 
effects as those required to match the observations requires the appearance in the dynamics of a strong turbulent regime, 
restated as a laminar one in the presence of a significant shear viscosity. Since this picture is not yet settled down (see |Balbus| 



2011), it is very important that an ideal hydrodynamical scheme on a Schwarzschild background like the one offered by the 



Polish doughnut is able to account for a material inf ailing onto a black hole. 

In this work we revisited the Polish doughnut model of accretion disks for a perfect fluid circularly orbiting around a 
Schwarzschild black hole with the effective potential approach for the exact gravitational and centrifugal contributions. We 
take advantage of the formal analogy between the fluid when the pressure vanishes and the test particle orbiting in the same 
background to get a comparison between the Polish doughnut, which is supported by the pressure, and the geodetic disk. Our 
analysis provides a revisited theoretical framework to characterize the accretion processes in presence of general relativistic 
effects. Indeed we formulate the Polish doughnut model in such a way that the fluid dynamics can be interpreted in terms of 
the fundamental stability properties of the circular orbits in Schwarzschild background. 

We extensively analyzed the Polish doughnut configurations for the fluid I and the particle L angular momenta, taken 
respectively constant throughout the entire toroidal surface. Then we propose a reinterpret at ion of torus physics, with respect 
to the its shape and equilibrium dynamics in terms of the parameter I (and L), and in terms of the parameter K that was 
naturally established by introducing the effective potential for the fluid motion. This new parameter K has been derived by 
exploiting the methodological and formal analogy with the effective potential approach to the test particles motion. Note that 
(/, K) fully describe the toroidal fluid configuration. This procedure made it possible to emphasize the pressure influence in 
the equilibrium dynamics with respect to the case of dust, the latter being treated as a set of test particles not subjected 
to pressure. This dual aspect, methodological and procedural, we thought it required a complete and deep analysis of the 
behavior of the momenta / and L on the plans and disk orbits, and of L as a function of /, and of the surface characteristic 
parameters as a function of K. 
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Figure 25. Upper panel: r± as a function of £. Gray region corresponds to Aq < 0, white region to Aq > 0. Aq = in r = r±. Inset: 
curves Aq = constant. The regions of positive and negative Aq are underlined. Lower panel: table summarizing the regions Aq > 
(light-gray) and A n < (gray). 

The main steps of this analysis are: 

Radial pressure gradient vs angular momentum C We studied the pressure radial gradient as a function of r and of the an- 
gular momentum C. The pressure is a decreasing function of r in r < ri co . and for C < Ck it is p < 0, this means that Ck 
identifies the pressure minimum points located in [n co ,n sco ], and pressure maximum points in r > ri sco (see Fig. [I] right 
panel). However, at fixed orbit r, Ck is always a minimum point of r ie at fixed r the pressure decreases until the angular 
momentum reaches the values C = Ck, and then increases with C (note that ri sco is a minimum point of Ck)- 



Angular momentum C vs. fluid angular momentum £ We found that at fixed £ the angular momentum C(r, I) has a maxi- 
mum for r = T- G [rico, risco], and a minimum in r+ > n sco , and C(r) increases for r < r_ and r > r+ (Fig. [3|. We identify 
three possibilities: £ < 3^/3/2, where the momentum C(r) increases with r; 3^3/2 < £ < 3\/3, where C increases with r up 
to the maximum point r_, decreases up to the minimum r+; £ > 3a/3, where C increases with r. 



Pressure gradients vs fluid angular momentum £ The fluid pressure decreases for r < rf and r > r+ as well as for £ < 3y3/2 
(Fig. [4]). The situation is much more complicated for fluid with higher angular momenta (£ > 3^/3/2). It appears necessary to 
consider fluids with 3^3/2 < £ < 3\/3 and with £ > 3\/3 separately. In the first case we observe the presence of a ring [r_, r+] 
where the fluid pressure increases with the orbital radius, r_ being a pressure minimum and r+ a pressure maximum. In the 
case £ > 3V3 we find two rings, [r l ,r_] and [r z + ,r+], where the fluid pressure is an increasing function of r. Figure [4] upper 
right describes this situation from a different point of view: for what concerns the variation of its hydrodynamic pressure with 
the orbits and the angular momentum, the fluid dynamics is basically split into two zones, r < ri co and r > n co , respectively: 
in the first region, the fluid pressure decreases with increasing angular momentum up to £ — £ r (minimum), then it increases 
with £ up to £k (maximum), finally it decreases with £. The trend is precisely the same in the region r > n co , but £k and £ r 
are now the minimum and maxima momentum respectively, i.e. the fluid pressure decreases with £ < £k, grows in the range 
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[IkAA and then decreases with £. The angular pressure gradient has been studied and compared to the radial and angular 
gradient in Sec. [4] 

The Boyer surfaces and polytropic equation of state In Sec. [5] we drew a complete and analytic description of the toroidal 
surface of the disk. A key role of this analysis was played by the effective potential approach: the toroidal disk can be described 
once one gives the effective potential and the fluid angular momentum. Our results concerning the disk shape and structure 
can be summarized as follows: 



(i) the distance from the source of the torus inner surface (S — y2 — 2M), increases with increasing angular momentum of 
the fluid but decreases with increasing energy function defined as the value of the effective potential for that momentum; 

(ii) the surface maximum height (torus thickness - h), increases with the energy and decreases with the angular momen- 
tum:the torus becomes thinner for high angular momenta, and thicker for high energies; 

(iii) the location of maximum thickness of the torus moves towards the external regions with increasing angular momentum 
and energy, until it reaches a maximum an then decreases; 

(iv) the surface maximum diameter (A) increases with the energy, but decreases with the fluid angular momentum; 

(v) the distance of the torus inner surface from the structure inner surface (5 = yi — yi) increases with the angular 
momentum and decreases with the energy; 

(vi) the distance of the structure inner surface from the horizon (5 = y± — 2M) increases with the energy and decreases 
with the fluid angular momentum. 



The accreting fluids with a polytropic equation of state were studied in the Section [6] divided into two classes identified by 
their polytropic index. 



The fluid angular velocity In Section [7] we analyzed the fluid proper angular velocity $ and we have compared the pr oper 
velocity with the Keplerian one &k in Sectionjs] the velocity $ < &k in r > r+ and for all r when £ < 3^/3/2 (see Fig. 22 ). As 
for the analysis of the angular momentum L(r, I), a distinction is made between fluids with £ £ [3 y/3/2, 3\/3], characterized by 
a ring [r_, r+] where the proper velocity is higher then the kepler one, and fluids with £ > 3\/3, where the proper fluid velocity 
is greater than the Keplerian one in the inner regions, r < r_. Finally, we have investigated the regions and momenta where 
the difference $ — &k is maximum (see Fig. |22[). We concluded the study of velocity fields by analyzing the fluid relativistic 
velocity Q and the Von Zeipel surfaces in Section [8] The velocity Q has a maximum at n co and increases with the increasing 
of £. We also studied the differences Aq = (Q — Q^)Ma between the relativistic fluid velocity and the Keplerian one Q K (see 
Fig. 23). At fixed angular momentum £ > 3^3/2, this difference has a maximum and a minimum, while for £ < 3^/3/2 it is 
always decreasing. The fluid velocity Q is lower then the Keplerian for r < r_ and r > r+, and for any r when £ < 3^/3/2, 
instead for fluid with with angular momentum greater then £ = 3^/3/2 there is a ring [r_,r+] where the fluid velocity Q is 



larger then the Keplerian one (see Fig. 25). 



The fundamental merit of the present work is that the analysis of the Polish doughnut features maintain the explicit 
presence of the effective potential in all the basic expressions describing the matter distribution. In fact, this allows to keep 
in direct contact our study with the behavior of free test particles moving the gravitational field of the central object and 
following circular orbit (the dust, pressure- free limit of the present analysis). What is significant here is the possibility to 
compare configuration of the considered fluid, as described by certain values of the parameters K and /, with the behavior of 
the test particle system, characterized by the same values of the corresponding parameters K and L. These latter quantities 
have a precise meaning for the particle (energy and angular momentum as viewed at infinite distance) , while the corresponding 
parameters for the fluid must, on this level, regarded as a classification criterion for the torus morphology. Thus retaining 
the effective potential in the Polish doughnut treatment allows to identify the Boyer surfaces in terms of parameters that 
have a precise meaning for a different, but closely related, context, the pressure-free fluid, i.e. the fundamental features of 
the Schwarzschild spacetime. If this study can not yet directly offer a paradigm for the comparison with the observations, 
nonetheless, it makes a concrete step in this direction. We are now able to constrain the morphology of the equilibrium 
configuration by fixing the two parameters K and /, i.e. specifying their value, we can define the basic features of the torus 
shape and of the velocity field in different space regions (this is synthetically sketched in the ten final remarks, achieved 
by our systematic analysis) and the comparison with the corresponding profile of the particle motion, where the physical 
comprehension is settled down, offer a valuable tool to interpret what the observations trace out. Our model will be completed 
including other relevant ingredients, like dissipative effects and the magnetic field and this is probably the natural development 
of the conceptual paradigm we fixed here. 
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Figure Al. (Color online) Left panel: conditions for G r > 0: ^ a < cr^ when < L < 2\/3M and r > ri co , and when L > 2\f?>M 
and ri co < r < r~ (a = 1) and r > r+(cr = 1). G r > for all < a < 1 when L > 2v / 3M and r~ (<r = 1) < r < r+(cr = 1), while it is 
always positive when r < r\ co . Right panel: conditions for G r < 0: cr^ < a < 1 when < L < 2\fZM and r > ri co , and when L > 2\f?>M 
and ri co < r < r^(a = 1) and r > r^(cr = 1). G r < for all < a < 1 when r < n co , while it is always negative when L > 2v / 3M and 
rj~(<7 = 1) < r < rj(cr = 1). Dotted lines mark r = ri co , dashed lines r = n sco , dot-dashed lines L = 2\f?>M . 



Table Al. Regions where d r L = 0. 
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APPENDIX A: THE RADIAL GRADIENT G R AS FUNCTION OF a. 



In Section [3.1 1 and |3.3| we detailed the study of the radial gradient G r as a function of the angular momentum L and of the 
fluid angular momentum Z with a re-parametrization that is independent from the equatorial plane cr. However, an explicit 
study of G r as a function of cr is useful for two main reasons: first, to have a direct comparison when the re-parametrization 
is not possible (as for the polar gradient G&). Secondly, the angular dependence is necessary to build a three-dimensional 
characterization of the torus. 

In Section |3] we showed that the critical points of the pressure are denned by the condition Eq. [14] In terms of the polar 
coordinate, it reads: 



L 

M 



M{r - 3M) 



(Al) 



The condition of existence of cr^ are r > n co for < L ^ 2y / 3M, and n co < r ^ (cr = 1) and r ^ rf (cr = 1) for L > 2y / 3M, 
where (cr — 1) are the circular orbit radii for a test particle in Schwarzschild spacetime denned by Eq. M5] evaluated at 
(7 = 1. The point ri sco is a maximum for cr^: the condition r£ = r\ — n sco is fulfilled when = a^* = \fl? /2\fZM . When 
2^M, it is ctk* 



Al 



where rf~ 



are introduced in Eq. 



L = 2V3M, it is = 1. The sign of G r is summarized in Figs 

We characterize now the angular momentum L as an explicit function of 
not denned in the interval [r[ 



19 



where ai r = yj 1 {r r3 — -. ai r is denned for r > 2M when < Z 

when Z > 3a/3M. r = r\ co is a maximum for cr^ r , where a/ r . = g\ 
The critical points of L are determined by the condition: 



and a. We know from Section [3. 2| that L is 
In terms of a this region corresponds to ^ a ^ crj r , 

< 3^3M, and for 2M < r < rf(cr = 1) and r > r+(cr = 1) 



3^3M 



When Z 3^3M, it is erf 1 



F( 



2M) 2 



Mr 3 



20). <7 K is denned for r > ri co when < I ^ 3y3/2M, and for ri co < r ^ r_(cr — 1) and r ^ r+(a 



< 3V3M, and also for r ^ r+(cr = 1) for , 



(see Eq. 

at a = 1. (Jk is maximum for r — n sco , where cr^ — cr : 

Finally the study of G r as a function of (r, Z, cr) is illustrated in Figs. A2 



^ 3^3M, where r±(cr = 1) are the radii denned by E q. 
= y/2a* . These results are summarized in Table 
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Figure A2. (Color online) Left panel: conditions for G r > 0: o\ r < a < cr^ when < I < 3y / 3/2M and r > n co , and when 
3^/3j2M < I < 3V3M and r ico < r < r_(cr = 1) and r > r + (cr = 1), and when I > 3\/3M and r > r + (cr = 1); a tr < a < 1 
when 3y/3/2M < Z < 3^3M and r_(cr = 1) < r <+ (cr = 1), and when I > 3^3M and rf(a = 1) < r < r+(cr = 1); a Y K < a < 1 
when Z > 3V3M and r~(cr = 1) < r < r_(cr = 1); a Y K < a < <j lr when I < 3V3M and 2M < r < r Zco , and when I > 3a/3M and 
2M < r < t z _ (<j = 1). Right panel: conditions for G r > 0: < a < ai r and crj^ < a < 1 when < Z < 3y / 3/2M and r > n co , and 
when 3^3/2M < Z < 3\/3M and r ico < r < r_(cr = 1) and r > r + (cr = 1), and when I > Sv^M and r > r + (cr = 1); ) < a < a tr 
when 3y / 3/2M < Z < Sv^M and r_(cr = 1) < r <+ (cr = 1), and when I > 3^3M and rf(a = 1) < r < r + (cr = 1); < a < a Y K when 
I > 3V3M and r~ (a = 1) < r < r_(cr = 1); < a < cr^ and cr Zr < cr < 1 when Z < 3V3M and 2M < r < n co , and when I > 3\/3M 
and 2M < r < r~ (cr = 1). G r < for all < a < 1 when L > 3\/3M and r_(cr = 1) < r < r+(cr = 1). Dotted lines mark I = 3\/3M, 



dot-dashed lines Z = 3y 3/2M, dashed lines r 
region Z/M G [6, 11] and r/M G [2,r Zco ]. 



risco and dotted gray lines r 
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Black regions are forbidden. Insets: zoom of the 



APPENDIX B: ANALYSIS OF THE PROPER VELOCITY PROFILE 

We redefine the radius in Eq. [97] in terms of a and Z: 



r x /M-- 
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for a / (for a = it is r J = 2M). 

The fluid proper angular velocity increases with r/M (<£' > 0) when: 
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Figure Bl. Left panel: conditions for < 0. Right panel: conditions for > 0. Solid line marks erf as a function of l/M, dashed line 
l/M = 3VS. 



(see Fig. Bl ). It follows that the angular velocity increases in 2M < r < r x < 3M until it reaches the point r = r x , that is a 
maximum, where 2.5M < r x < 3M with angular momentum l x on the equatorial plane <j x , then it decreases with the radius 



(see Fig. Bl). 
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